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ABSTRACT 

Aims. We present an asteroseismological study on the two high-gravity pulsating PGl 159 (GW Vir or DOV) stars, PG 2131-1-066 and 
PG 1707H-427, and on the pulsating [WCE] star NGC 1501. All of these stars have been intensively scrutinized through multi-site 
observations, so they have well resolved pulsation spectra. 

Methods. We compute adiabatic g-mode pulsation periods on PGl 159 evolutionary models with stellar masses ranging from 0.530 to 
O.74IM0. These models take into account the complete evolution of progenitor stars, through the thermally pulsing AGB phase, and 
born-again episode. We constrain the stellar mass of PG 2131+066, PG 1707-1-427, and NGC 1501 by comparing the observed period 
spacing with the asymptotic period spacing and with the average of the computed period spacings. We also employ the individual 
observed periods in search of representative seismological models for each star. 

Results. We derive a stellar mass of 0.627 Mq for PG 2131-t-066, 0.597 Mq for PG 1707-1-427, and 0.571 Mq for NGC 1501 from 
a comparison between the observed period spacings and the computed asymptotic period spacings, and a stellar mass of 0.578 Mq 
for PG 21314-066, 0.566 Mg for PG 1707+427, and 0.576 Mq for NGC 1501 by comparing the observed period spacings with the 
average of the computed period spacings. We also find, on the basis of a period-fit procedure, asteroseismological models represen- 
tatives of PG 2131+066 and PG 1707+427. These best-fit models are able to reproduce the observed period patterns of these stars 
with an average of the period difi^erences of (511, = 1.57 s and (511, = 1.75 s, respectively. The best-fit model for PG 2131+066 has an 
effective temperature T^ff = 102 100 K, a stellar mass M, = 0.589 Mq, a surface gravity logg = 7.63, a stellar luminosity and radius 
of log(L»/Lo) = 1.57 and \og(R,/R0) = -1.71, respectively, and a He-rich envelope thickness of Me„v = 1-6 x 10"^ Mq. We derive a 
seismic distance d ~ 830 pc and a parallax n ~ 1.2 mas. The best-fit model for PG 1707+427, on the other hand, has T^g = 89 500 
K, M, = 0.542 Mg, logg = 7.53, log(L./Lo) = 1.40, \og(R,/RQ) = -1.68, and M^„^ = 2.5 x 10"^ Mq, and the seismic distance and 
parallax are d ~ 730 pc and n ~ 1.4 mas. Finally, we have been unable to find an unambiguous best-fit model for NGC 1501 on the 
basis of a period-fit procedure. 

Conclusions. This work closes our short series of asteroseismological studies on pulsating pre-white dwarf stars. Our results demon- 
strate the usefulness of asteroseismology for probing the internal structure and evolutionary status of pre-white dwarf stars. In par- 
ticular, asteroseismology is able to determine stellar masses of PGl 159 stars with an accuracy comparable or even better than spec- 
troscopy. 
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1. Introduction 

Pulsating PGl 159 stars (also called GW Vir or DOV stars) 
are very hot hydrogen-deficient post- Asymptotic Giant Branch 
(AGB) stars with surface layers rich in helium, carbon, and oxy- 
gen (Werner & Herwig 2006) that exhibit multiperiodic lumi- 
nosity variations with periods ranging from 5 to 50 minutes, at- 
tributable to non-radial pulsation g-modes (see Winget & Kepler 
2008 for a recent review). PGl 159 stars are thought to be the 
evolutionary link between Wolf-Rayet type central stars of plan- 
etary nebulae and most of the hydrogen-deficient white dwarfs 
(Althaus et al. 2005). It is generally accepted that these stars 
have their origin in a born-again episode induced by a post- AGB 
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helium thermal pulse — see Iben et al. (1983), Herwig et al. 
(1999), Lawlor & MacDonald (2003), and Althaus et al. (2005) 
for recent references. 

Recently, considerable observational eff'ort has been invested 
to study pulsating PGl 159 stars. Particularly noteworthy are 
the works of Vauclair et al. (2002) on RX J21 17.1+3412, Fu 
et al. (2007) on PG 0122+200, and Costa et al. (2008) and 
Costa & Kepler (2008) on PG 1 159-035. These stars have been 
monitored through long-term observations carried out with the 
Whole Earth Telescope (Nather et al. 1990). On the theoret- 
ical front, recent important progress in the numerical model- 
ing of PGl 159 stai-s (Althaus et al. 2005; Miller Bertolami & 
Althaus 2006, 2007a) has paved the way for unprecedented as- 
teroseismological inferences for the mentioned stars (Corsico & 
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Althaus 2006, Corsico et al. 2007a, 2007b, 2008). The new gen- 
eration of PG1159 evolutionary models of Miller Bertolami & 
Althaus (2006) is derived from the complete evolutionary his- 
tory of progenitor stars with different stellar masses and an elab- 
orate treatment of the mixing and extra-mixing processes during 
the core helium burning and born-again phases. The success of 
these models at explaining the spread in surface chemical com- 
position observed in PGl 159 stars (Miller Bertolami & Althaus 
2006), the short born-again times of 4334 (Miller Bertolami & 
Althaus 2007b), and the location of the GW Vir instability strip 
in the log T^ff - \ogg plane (Corsico et al. 2006) renders relia- 
bility to the inferences drawn from individual pulsating PGl 159 
stars. 

Besides the mentioned three well-studied pulsating PGl 159 
stars, there exist two other variable stars of this class that have 
been also intensively scrutinized through the multi-site obser- 
vations of the WET: PG 2131+066 and PG 1707+427. In addi- 
tion, there is a variable central star of planetary nebula (PNNV), 
NGC 1501, which has been the subject of a nearly continu- 
ous photometric coverage from a global observing campaign by 
Bond et al. (1996). We briefly summarize the properties of these 
stars below. 

PG 2131+066 was discovered as a variable star by Bond et 
al. (1984) with periods of about 414 and 386 s, along with some 
other periodicities. On the basis of an augmented set of peri- 
ods from WET data, Kawaler et al. (1995) obtained a precise 
mass determination of M» - 0.61 Mq, a luminosity of 10 L©, 
and a seismological distance from the Earth of c/ = 470 pc. 
Spectroscopic constraints of Dreizler & Heber (1998), on the 
other hand, gave M, = 0.55 Mq, T^ff = 95 000 K, 39.8 Lq, and 
log^ = 7.5 for PG 2131+066. By using this updated determina- 
tion of the effective temperature. Reed et al. (2000) refined the 
procedure of Kawaler et al. (1995) and found M» = 0.608 Mq, 
L, - 26 Lq and d = 668 pc. 

PG 1707+427 was discovered to be a pulsator by Bond et 
al. (1984). Dreizler & Heber (1998) obtained T^g = 85 000 K, 
logg = 7.5, and then M, = 0.54 Mq and L, = 25 Lq were 
inferred from their spectroscopic study. Recently, Kawaler et 
al. (2004) reported the analysis of multi-site observations of 
PG 1707+427 obtained with WET. Preliminary seismic analysis 
by using 7 independent { - I modes with periods between 334 
and 910 s suggest an asteroseismological mass and luminosity 
of 0.57 Mq and 23 L©, respectively. 

NGC 1501 was classified as a [WCE] star, an early low-mass 
Wolf Rayet-type PNNV with spectra dominated by strong he- 
lium and carbon emission lines (Werner & Herwig 2006). The 
effective temperature and gravity of this star are Teff = 1 34 000 
K and log^ = 6.0 (Werner & Herwig 2006). The variable nature 
of NGC 1501 was discovered by Bond & Ciardullo (1993). The 
star shows ten periodicities ranging from 5200 s down to 1154 
s, although the largest amplitude pulsations occur between 1 154 
s and 2000 s. Based on period-spacing data. Bond et al. (1996) 
found a stellar mass of 0.53 + 0.03 Mq for NGC 1501. 

In this work we complete our small survey of astero- 
seismological inferences on pulsating PGl 159 stars — see 
Corsico et al. (2007a, 2007b, 2008) for the previous studies 
of this series — by performing a detailed study of the GW 
Vir stars PG 2131+066 and PG 1707+427, and the [WCE] star 
NGC 1501. We employ the same stellar models and numeri- 
cal tools as in our previous works. In particular, we go beyond 
the mere use of the observed period-spacing data by perform- 
ing, in addition, detailed period-to-period fits on the pulsation 
spectrum of these stars. In our approach, we take full advan- 
tage of the state-of-the-art PGl 159 evolutionary models devel- 



oped by Miller Bertolami & Althaus (2006). The paper is orga- 
nized as follows. In Section |2] we briefly describe our PGl 159 
evolutionary models. In Sect. [3] we derive the stellar mass of 
PG 2131+066, PG 1707+427, and NGC 1501 by using the ob- 
served period-spacing data alone. In Sect. |4] we infer structural 
parameters of these stars by employing the individual observed 
periods. In this section we derive asteroseismological models 
representative of PG 2131+066 and PG 1707+427 dO), and 
discuss their main structural and pulsational characteristics ( I4.2l i. 
In Sect. |5] we compare the results of the present paper with those 
of the asteroseismological study of Corsico & Althaus (2006) 
(hereinafter CA06). Finally, in Sect. |6] we summarize our main 
results and make some concluding remarks. 

2. Evolutionary models and numerical tools 

The pulsation analysis presented in this work relies on a new 
generation of stellar models that take into account the complete 
evolution of PGl 159 progenitor stars. Specifically, the stellar 
models were extracted from the evolutionary calculations re- 
cently presented by Althaus et al. (2005), Miller Bertolami & 
Althaus (2006), and Corsico et al. (2006), who computed the 
complete evolution of model star sequences with initial masses 
on the ZAMS ranging from 1 to 3.75 Mq. All of the post-AGB 
evolutionary sequences were computed using the LPCODE evo- 
lutionary code (Althaus et al. 2005) and were followed through 
the very late thermal pulse (VLTP) and the resulting born-again 
episode that gives rise to the H-deficient, He-, C-, and 0-rich 
composition characteristic of PGl 159 stars. The masses of the 
resulting remnants are 0.530, 0.542, 0.556, 0.565, 0.589, 0.609, 
0.664, and 0.741 Mq. For details about the input physics and evo- 
lutionary code used, and the numerical simulations performed to 
obtain the PGl 159 evolutionary sequences employed here, we 
refer the interested reader to the works by Althaus et al. (2005) 
and Miller Bertolami & Althaus (2006, 2007a). 

We computed f = 1 g-mode adiabatic pulsation periods with 
the same numerical code and methods we employed in our pre- 
vious works, see Corsico & Althaus (2006) for details. In ad- 
dition, we performed nonadiabatic computations with the help 
of the code employed in Corsico et al. (2006) to evaluate the 
pulsational stability of the asteroseismological models presented 
in Sect. H] We analyzed about 3000 PGl 159 models covering 
a wide range of effective temperatures (5.4 ^ log(reff) ^ 4.8) 
and luminosities (0 ^ log{Lt/LQ) ^ 4.2), and a range of stellar 
masses (0.530 < MJMq < 0.741). 

3. Mass determination from the observed period 
spacing 

In this section we constrain the stellar mass of PG 2131+066, 
PG 1707+427, and NGC 1501 by comparing the asymptotic pe- 
riod spacing and the average of the computed period spacings 
with the observed period spacing. These approaches take full ad- 
vantage of the fact that the period spacing of PGl 159 pulsators 
depends primarily on the stellar mass, and the dependence on 
the luminosity and the He-rich envelope mass fraction is neg- 
Hgible (Kawaler & Bradley 1994; Corsico & Althaus 2006). 
Most of the published asteroseismological studies on PGl 159 
stars rely on the asymptotic period spacing to infer the total 
mass of GW Vir pulsators, the notable exception being the works 
by Reed et al. (2000) for PG 2131+066, Corsico et al. (2007a) 
for RX J21 17.1 +3412, Corsico et al. (2007b) for PG 0122+200, 
and Kawaler & Bradley (1994) and Corsico et al. (2008) for 
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Table 1. Stellar masses for all of the intensively studied pulsating PG 1159 stars, including also one pulsating [WCE] star All 
masses are in solar units. 



Star 


An? 




Approximate 
formula 


Period fit 


Pulsations 
(other works) 


Spectroscopy 


NGC 1501 


0.571" 


0.576" 


0.530" 




0.55J (asymptotic analysis) 


0.56 


RX J2117.1+3412 


0.5681^ 


0.560'' 


0.525" 


0.565'' 


0.56'' (asymptotic analysis) 


0.72 


PG 1159-035 


0.577-0.585'' 


0.56^ 


0.570" 


0.565'' 


0.59' (asymptotic analysis) 


0.54 


PG 2131+066 


0.627" 


0.578" 


0.609" 


0.589" 


0.61"= (period fit) 


0.55 


PG 1707+427 


0.597" 


0.566" 


0.587" 


0.542" 


0.57° (asymptotic analysis) 


0.53 


PG 0122+200 


0.625*^ 


0.567^ 


0.593" 


0.566"^ 


0.69*^ (asymptotic analysis) 


0.53 



Notes: "This work. ''Corsico et al. (2007a). '^Corsico et al. (2007b). ''Corsico et al. (2008). "^Reed et al. (2000). 'Fu et al. (2007). ^Kawaler et al. 
(2004). ''Vauclair et al. (2002). 'Costa et al. (2008). JBond et al. (1996). 
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Fig. 1. The dipole asymptotic period spacing in terms of the ef- 
fective temperature. Numbers along each curve denote the stellar 
mass (in solar units). Dashed (solid) lines correspond to evolu- 
tionary stages before (after) the turning point at the maxima ef- 
fective temperature of each track. Also shown are the locations 
of PG 2131+066, PG 1707+427, and NGC 1501. The masses 
of these stars as derived by comparing AIl^^j with AII^^j are 
M. = 0.627 Mq, M. = 0.597 Mo, and M. = 0.571 Mq, respec- 
tively. Note that to infer the mass of NGC 1501, we have consid- 
ered the stages before the evolutionary knee (dashed lines). 



et al. 1996). The asymptotic period spacing is computed as 
An^ = Ho/ ^/I(JTT), where 



f 



iN/r)dr 



(1) 



and is the Brunt-Vaisala frequency (Tassoul et al. 1990). 
From the comparison between the observed An°^ J andAH^^j we 
found a stellar mass of 0.627 Mq for PG 2131+066, 0.597 Mq 
for PG 1707+427, and 0.571 Mq for NGC 1501 (second column 
in Table [B- 

The method employed here is computationally inexpensive 
and widely used because it does not involve pulsational cal- 
culations. However, we must emphasize that the derivation of 
the stellar mass using the asymptotic period spacing may not 
be entirely reliable in pulsating PG1159 stars that pulsate with 
modes characterized by low and intermediate radial orders (see 
Althaus et al. 2007). This is particularly true for PG 2131+066 
and PG 1707+427. This shortcoming of the method is due to 
that the asymptotic predictions are strictly valid in the limit of 
very high radial order (long periods) and for chemically homo- 
geneous stellar models, while PG1159 stars are supposed to be 
chemically stratified and characterized by strong chemical gra- 
dients built up during the progenitor star life. A more realistic 
approach to infer the stellar mass of PG1159 stars is presented 
below. 



PG 1159-035. To assess the total mass of NGC 1501 we have 
considered the high-luminosity regime of the evolutionary se- 
quences, while forPG 2131+066 and PG 1707+427 we have fo- 
cused on the stages following the "evolutionary knee" for the PG 
1 159 stars, i.e. the low-luminosity regime. 



3.1. First method: comparing fhe observed period spacing 
(AIl'^) witii tine asymptotic period spacing fAII^j 

Fig. [T] displays the asymptotic period spacing for {. - \ modes 
as a function of the effective temperature for different stel- 
lar masses. Also shown in this diagram is the location of 
PG 2131+066, with T^s = 95 + 5 kK (Dreizler & Heber 1998), 
and An° J = 21 .6 + 0.4 s (Reed et al. 2000), PG 1707+427, with 
Teff = 85+5 kK (Dreizler & Heber 1998), and An°, = 23.0+0.3 
s (Kawaler et al. 2004), and NGC 1501, with T^s = 134 + 5 
kK (Werner & Herwig 2006), and An° j = 22.3 + 0.3 s (Bond 



3.2. Second metliod: comparing fhe observed period 

spacing (AYl'^) witli tlie average of tiie computed period 

spacings (AUc) 

The average of the computed period spacings is assessed as 
Allf - {N - 1)"' 2a Alli, where the "forward" period spacing 
is defined as Alli = ^k+\ - (k being the radial order) and 
is the number of computed periods laying in the range of the 
observed periods. ForPG 2131+066, U,, e [340,600] s, accord- 
ing to Kawaler et al. (1995); for PG 1707+427, H^ e [330, 920] 
s, according to Kawaler et al. (2004); and for NGC 1501, Uk e 
[1 150, 2000] s, according to Bond et al. (1996). 

This method is more reliable for the estimation of the stellar 
mass of PG1159 stars than that described above because, pro- 
vided that the average of the computed period spacings is eval- 
uated at the appropriate range of periods, the approach is ap- 
propriate for the regimes of short, intermediate and long peri- 
ods (i.e., VA;) as well. When the average of the computed period 
spacings is taken over a range of periods characterized by high k 
values, then the predictions of the present method become closer 
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Fig. 2. Same as Fig. [T] but for the average of the computed pe- 
riod spacings. For PG 2131+066 and PG 1707+427 only the 
stages after the "evolutionary knee" have been plotted. The 
masses of PG 2131+ 066, PG 1707+427, and NGC 1501 as de- 
rived by comparing Allf^i with An°^j are M, - 0.578 Mq, 
M« = 0.566 Mq, and M» = 0.576 Mq, respectively. As in Fig. 
[U we infer the mass of NGC 1501 by considering the stages be- 
fore the evolutionary knee (dashed lines). 



to those of the asymptotic period spacing approach. On the other 
hand, the present method requires of detailed period computa- 
tions, at variance with the method described in the above sec- 
tion. In addition, we note that both methods for assessing the 
stellar mass rely on the spectroscopic effective temperature, and 
the results are unavoidably affected by its associated uncertainty. 

In Fig. |2] we show the run of average of the computed pe- 
riod spacings ({ = 1) for PG 2131+066, PG 1707+427, and 
NGC 1501 in terms of the effective temperature for all of our 
PGl 159 evolutionary sequences. The run of Allf depends on the 
range of periods on which the average of the computed period 
spacing is done. Note that the lines shown in Fig. |2] are very 
jagged and jumped. This is because that, for a given star, the av- 



erage of the computed period spacings is evaluated for a fixed 
period interval, and not for a fixed fc-intervalQ. 

By adopting the effective temperature of PG 2131 +066, 
PG 1707+427, and NGC 1501 as given by spectroscopy we 
found a stellar mass of 0.578 Mq, 0.566 Mq, and 0.576 Mq, re- 
spectively. Our results are shown in the third column of Table 
[U These values are 8.5% (for PG 2131+066) and 5.5% (for 
PG 1707+427) smaller than those derived through the asymp- 
totic period spacing, showing once again that the asymptotic ap- 
proach overestimates the stellar mass of PGl 159 stars that, like 
PG 2131+066 and PG 1707+427, exhibit short and intermedi- 
ate pulsation periods (see Althaus et al. 2008a). On the con- 
trary, there is a very small discrepancy (in the opposite direc- 
tion) of ~ 0.9% for the case of NGC 1501, showing that in the 
long-period regime the results for the stellar mass obtained using 
the asymptotic period spacing and the average of the computed 
period spacings nicely agree each other. A similar situation is 
found in the case of RX J21 17. 1 +3412 (Corsico et al. 2007a). 

3.3. Third method: using an approximate formula 

To compare with previous works, we make an additional es- 
timation of the stellar mass of PG 2131+066, PG 1707+427, 
and NGC 1501 using the approximate expression for the overall 
structure parameter Hq derived by Kawaler & Bradley (1994): 



Ho ~ 15.5 



Mr 



\100Lf 



-0.035 



-0.00012 



(2) 



where qy is the He-rich envelope mass fraction. This expression 
is derived by considering the dependence of the asymptotic pe- 
riod spacing on the total mass, stellar luminosity, and thickness 
of the He-rich outer envelope for a large grid of "quasi evolution- 
ary" PGl 159 models in the luminosity range 1.6 ^ logiLt/Lo) :$ 
3.0. Both the present method and the method described in Sect. 
13. H are almost equivalent because they are based on the asymp- 
totic period spacing. 

Due to the very weak dependence of Ho on qy, we arbitrar- 
ily fix it to a value of 10"^. Since the luminosity is not known 
at the outset, we can compute it as L, = AncrRlT*^, where 
RI - GMt/g. We use the values of g and Tetf inferred through 
spectroscopy. Assuming that Do is known from the observed pe- 
riod spacing (Ho ~ V^(^ + 1) An°), we obtain an estimation 
of the stellar mass from Eq. (|2]l. Our results are shown in the 
fourth column of Table[I] For PG 2131+066 and PG 1707+427 
the stellar masses obtained in this way are in very good agree- 
ment with our values derived from the asymptotic period spacing 
(first row in Table [T]i. This is not an unexpected result because, 
as mentioned, the expression of Kawaler & Bradley (1994) is 
also based on the asymptotic period spacing. The slight diff'er- 
ences found (below » 2.5%) could be due to differences in the 
modeling of PG 1159 stars. For NGC 1501, instead, there is a 
substantial difference (» 7%) between the prediction of this for- 
mula and our value inferred from the asymptotic period spacing. 
This could be due to the inadequacy of the formula of Kawaler 



' As the star evolves towards higher (lower) effective temperatures, 
the periods generally decrease (increase) with time. At a given T^ff, 
there are A' computed periods laying in the chosen period interval. Later 
when the model has evolved enough (heated or cooled) it is possible that 
the accumulated period drift nearly matches the period separation be- 
tween adjacent modes (\Ak\ = 1). In these circumstances, the number of 
periods laying in the chosen (fixed) period interval is N ± 1, and Allf 
exhibits a little jump. 
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& Bradley (1994) for the high-luminosity regime characterizing 
the evolutionary status of NGC 1501. 

4. Constraints from the individual observed periods 

In this approach we seek pulsation models that best match the in- 
dividual pulsation periods of PG 2131+066, PG 1707+427, and 
NGC 1501. For the three stars, we assume that all of the ob- 
served periods correspond to ^ = 1 modes because the observed 
period spacings and the frequency splittings by rotation are con- 
sistent with £ - 1 (Kawaler et al. 1995; Kawaler et al. 2004; 
Bond et al. 1996). To measure the goodness of the match be- 
tween the theoretical pulsation periods (IlJ) and the observed in- 
dividual periods (11°), we follow the samex^ procedure as in our 
previous works — see, e.g, Corsico et al. (2007a) and Townsley 
et al. (2004). Specifically, we employ the quality function de- 
fined as 

1 " 

x\M,, Teff) = - 2 ™" [(n^ - n^)2] (3) 

/=1 

where n is the number of observed periods. The observed periods 
are shown in the first column of Tables|2]and[3]for PG 2131 +066 
and PG 1707+427, and in Table 7 of Bond et al. (1996) for 
NGC 1501. 

Next, we briefly explain the procedure we follow to found 
a model representative of a target star. For a given model (char- 
acterized by a given T^ff) corresponding to an evolutionary se- 
quence of stellar mass M,, we consider the first observed period 
of the list, namely n°, and compute the successive squared dif- 
ferences (n° - nj)^, where the radial order k varies from 1 to 
a given maximum value, k^n-^j^, which corresponds to a theoreti- 
cal period far longer that the maximum period observed in the 
star. Next, we retain the minor squared difference, and then we 
repeat the procedure but this time considering the second ob- 
served period, namely . After the minor squared difference 
associated with this observed period is stored, we proceed with 
the next observed period, and so until the minor squared differ- 
ence associated with the last observed period (11°) is stored. The 
next step is to calculate the sum of these differences and then 
obtain the value of x^(Mt, T^ff) for the model under considera- 
tion. It is worth mentioning that with this algorithm, the value 
of x^{Mt, Teff) does not depend on the particular order in which 
the observed periods are fitted. The complete algorithm is re- 
peated for all of the models of the sequence, and then, a curve 
of ;^'^(M,, Teff) versus Tgff is obtained for the complete sequence. 
This procedure is carried out for all of our sequences. For each 
star of interest, the PG 1 159 model that shows the lowest value 
ofx^ is adopted as the "best-fit model". 

4.1. The search for the best-fit models 

We evaluate the function x^ - X^i^*^ T^s) for stellar masses 
of 0.530, 0.542, 0.556, 0.565, 0.589, 0.609, 0.664, and 0.741Mo. 
For the efifective temperature we employed a much finer grid 
(Areft = 10-30 K). The quantity (^-y^ in terms of the effec- 
tive temperature for different stellar masses is shown in Fig. [3] 
for PG 2131+066 (upper panel), PG 1707+427 (middle panel), 
and NGC 1501 (lower panel), together with the correspond- 
ing spectroscopic effective temperatures. We prefer to show in 
our plots the quantity (x^Y"^ instead x^ in order to emphasize 
the location of models providing good agreements between ob- 
served and theoretical periods. As mentioned, the goodness of 
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Fig. 3. The inverse of the quality function of the period fit in 
terms of the effective temperature (see text for details). The ver- 
tical grey strip indicates the spectroscopic Teff and its uncertain- 
ties. The curves have been arbitrarily shifted upward (with a step 
of 0.1). Upper panel corresponds to PG 2131+066, middle panel 
corresponds to PG 1707+427, and lower panel corresponds to 
NGC 1501. For NGC 1501, the stages before the evolutionary 
knee are displayed with dotted lines. 



the match between the observed and theoretical periods is re- 
flected by the value of x^- The lower the value of x^, the bet- 
ter the period match. We will consider — admittedly somewhat 
arbitrarily — that a peak in the quality function with x^ (that 
is, (x^y^ > 0.2) is a good match between the theoretical and the 
observed periods. 

For PG 2131+066 we find one strong maximum of (x^Y^ 
for a model with M, = 0.589 Mq and T^s ^ 102 kK. Such 
a pronounced maximum in the inverse of x^ implies an ex- 
cellent agreement between the theoretical and observed peri- 
ods. Another much less pronounced maxima, albeit at effec- 
tive temperatures closer to the spectroscopic estimation for 
PG 2131+066, are encountered for M, = 0.609 Mq at T^s ~ 
92.4 kK and = 0.565 Mq at T^s ~ 99.8 kK. However, be- 
cause the agreement between observed and theoretical periods 
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for these models are substantially poorer than for the one with 
Mt - 0.589 Mq, we adopt this last model as the best-fit aster- 
oseismological model. A detailed comparison of the observed 
{ - l,m = periods in PG 2131+066 with the theoretical peri- 
ods of the best-fit model is provided in Table|2l The high quality 
of our period fit is quantitatively reflected by the average of the 
absolute period differences 



2 " 



(4) 



1=1 

where 611, = IT 



° - nT and by the root-mean-square residual 



(5) 



Note that cr^ is simply the function -\/x^ evaluated at the best-fit 

model. For the best-fit model of PG 2131+066 we obtain 611, - 
1.57 s and cr^j, = 2.32 s, which are indeed very small. 

For PG 1707+427 we find one strong maximum of (x^y^ 
for a model with M, = 0.542 Mq and T^ff ^ 89.5 kK, an eff^ec- 
tive temperature compatible with the spectroscopic determina- 
tion. Another somewhat less pronounced maximum is found for 
a model with M, = 0.556 Mq and ~ 83.6 kK. Despite the 
fact that this model has an effective temperature very close to the 
spectroscopic one, we choose the model with M» = 0.542 Mq 
as the best-fit model for PG 1707+427, because the period fit is 
characterized by a better quality. Table[3]shows a comparison be- 
tween observed ( = l,m - and computed periods of the best- 
fit model. We found in this case dYl, - 1.75 s and cr^„ = 1.99 
s. 

The situation for NGC 1501 is markedly different than for 
PG 2131+066 and PG 1707+427. Indeed, as clearly shown in 
the lower panel of Fig. [3] for this star the (x^y^ function ex- 
hibits numerous local maxima at several values of the effective 
temperature and the stellar mass (M»/Mq = 0.556, 0.565, 0.589, 
0.664, and 0.741) that have roughly the same amplitudes, mak- 
ing virtually impossible to isolate a clear and unambiguous seis- 
mological solution. Thus, for NGC 1501 we are unable to find 
a best-fit seismological model. This could be, in part, due to the 
fact that the periods detected in NGC 1501 can be associated to 
eigenmodes with radial orders k quite different from each other 
(with a mean spacing of Ak ^ 6), in such a way that it is easy to 
find numerous models (characterized by strongly different T^ff 
and Mt) that reproduce to a some extent the observed period 
spectrum of NGC 1501. It could also be that the impossibility 
to find a best fit model would be reflecting a different evolution- 
ary history for NGC 1501 than that assumed in this work for our 
PG1159 sequences. 

The fourth column in Tables |2] and [3] shows the rates of pe- 
riod change associated with the fitted modes for PG 2131+066 
and PG 1707+427, respectively. Our calculations predict that all 
of the pulsation periods increase with time (11 > 0), in accor- 
dance with the decrease of the Brunt- Vaisala frequency in the 
core of the models induced by cooling. At the effective tempera- 
tures of PG 2 1 3 1 +066 and PG 1 707 +427, cooling has the largest 
effect on 11, while gravitational contraction, which should result 
in a decrease of periods with time, becomes negligible and no 
longer affects the pulsation periods, except for the case of modes 
trapped in the envelope. Until now, the only secure measure- 
ments of n in pre-white dwarf stars are those of PG 1 159-035, 
the prototype of the class, by Costa et al. (1999) and more re- 
cently by Costa & Kepler (2008). In this last paper the authors 



Table 2. Observed and theoretical ({ — I) periods of the best- 
fit model for PG 2131+066 (M. = 0.589Mq, T^s = 102 179 
K, log(L»/LQ) = 1.57). Periods are in seconds and rates of pe- 
riod change (theoretical) are in units of 10"'^ s/s. The observed 
periods are taken from Kawaler et al. (1995). 







k 




unstable 


341.45 


341.88 


14 


2.04 


yes 




363.66 


15 


3.90 


yes 


384.27 


383.94 


16 


2.60 


yes 


403.93 


403.30 


17 


3.30 


yes 


426.36 


426.48 


18 


4.14 


yes 


450.28 


445.31 


19 


2.89 


yes 


462.39 


465.91 


20 


4.53 


yes 




488.81 


21 


4.11 


no 


507.91 


507.91 


22 


3.79 


no 



Table 3. Same as Table |2l but for the best-fit model for PG 
1707+427 (M. = 0.542 Mq, = 89 504 K, log(L./LQ) = 
1.40). The observed periods are taken from Kawaler et al. 
(2004). 



n° 




k 




unstah 


334.62 


333.26 


12 


2.13 


yes 




355.73 


13 


1.19 


yes 




378.53 


14 


2.17 


yes 




401.64 


15 


1.51 


yes 




423.68 


16 


1.84 


yes 


448.07 


447.31 


17 


2.46 


yes 




469.26 


18 


1.49 


yes 


494.39 


492.56 


19 


2.76 


yes 




516.50 


20 


2.61 


yes 


536.41 


537.67 


21 


1.66 


yes 




562.32 


22 


3.30 


yes 




584.78 


23 


2.37 


yes 




606.62 


24 


2.73 


yes 




632.11 


25 


3.25 


yes 




653.15 


26 


2.99 


yes 


677.89 


677.38 


27 


3.43 


yes 




701.20 


28 


3.03 


yes 


726.02 


722.82 


29 


3.95 


no 


745.78 


747.58 


30 


3.42 


no 




770.93 


31 


3.71 


no 




793.13 


32 


4.45 


no 




817.86 


33 


3.61 


no 




841.70 


34 


4.64 


no 




863.34 


35 


3.91 


no 




888.59 


36 


5.03 


no 


909.05 


912.34 


37 


4.21 


no 



obtained a mix of positive and negative values of H, indicating 
that for that star gravitational contraction is still important. In 
principle, the needed time interval that the observational data 
should cover in order to reach a measurement of a O in PG 1 159 
stars like PG 2131+066 and PG 1707+427 is of about ten years. 
Unfortunately, no future observations in the short term that could 
allow a determination of Et for these stars and thus to check the 
predictions of our models are foreseen. 

Finally, the last column in Tables |2] and [3] gives information 
about the pulsational stability/instability nature of the modes as- 
sociated with the periods fitted to the observed ones. Full nona- 
diabatic calculations employing the pulsation code described in 
Corsico et al. (2006) indicate that, for the case of PG 2 1 3 1 +066, 
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all the fitted modes except one (that with period 507.9 s) are pre- 
dicted to be unstable. For the case of PG 1707+427, our nona- 
diabatic computations are able to explain the existence of peri- 
odicities in the range 330 ^ 11° i 680 s only, while they fail to 
predict pulsational instability for the observed modes with peri- 
ods at 726, 746, and 909 s. 



4.2. Characteristics of the best-fit models for PG 2131+066 
and PG 1 707+427 

The main features of our best-fit model for PG 2131+066 are 
summarized in Table |4] where we also provide the parameters 
of the star extracted from other published studies. Specifically, 
the second column corresponds to spectroscopic results from 
Werner & Herwig (2006), whereas the third and fourth columns 
present results from the pulsation studies of Kawaler et al. 
(1995) and Reed et al. (2000), and from the asteroseismologi- 
cal model of this work, respectively. The number in parenthe- 
sis is the spectroscopic estimation of the stellar mass employing 
the evolutionary tracks of Miller Bertolami & Althaus (2006). 
In the present work, errors in T^s and logiLt/Lo) are estimated 
from the width of the maximum in the function with re- 
spect Teff and log(L,/L0), respectively. The error in the stellar 
mass comes from the grid resolution in M,. Errors in the rest 
of the quantities are derived from these values. The effective 
temperature of our best-fit model (Teff - 102 180 K) is some- 
what higher than — but still compatible with — the spectroscopic 
value (T^ff = 95 000±5 000 K). On the other hand, the total mass 
of the best-fit model (M* = 0.589 Mq) is in agreement with the 
value derived from the average of the computed period spac- 
ings (M. ~ 0.578 Mq), but at odds (~ 6% smaller) with that 
inferred from the asymptotic period spacing (M» = 0.627 Mq) 
(see Table[T]i. Also, the M» value of our best-fit model is substan- 
tially larger than the spectroscopic mass of 0.55 Mq derived by 
Miller Bertolami & Althaus (2006), but very similar to 0.58 Mq 
according to Werner & Herwig (2006). A discrepancy between 
the asteroseismological and the spectroscopic values of M» is 
generally encountered among PGl 159 pulsators — see Corsico 
et al. (2006, 2007 ab). Until now, the asteroseismological mass of 
PG 2 1 3 1 +066 has been about 1 1 % larger (AM. * 0.06 Mq) than 
the spectroscopic mass if we consider the early estimation for the 
seismological mass quoted by Reed et al. (2000) and the deriva- 
tion of Miller Bertolami & Althaus (2006) for the spectroscopic 
masfl In light of the best-fit model derived in this paper, this 
discrepancy is slightly reduced to about 7% (AM, ^ 0.04 Mq). 
Finally, our best-fit model for PG 2131+066 is somewhat more 
luminous and less compact than what is suggested by the results 
of Reed et al. (2000). 

The main properties of our best-fit model for PG 1707+427 
are shown in Table |5] The second column corresponds to spec- 
troscopic results from Werner & Herwig (2006), whereas the 
third and fourth columns present results from the pulsation study 
of Kawaler et al. (2004) and from the asteroseismological model 
of this work, respectively. As for the case of PG 2131+066, the 
effective temperature of our best-fit model for PG 1707+427 is 
slightly larger than the spectroscopic measurement, but even in 
good agreement with it. Regarding the stellar mass, our best- 



Table 4. The main characteristics of PG 2131+066. 



^ We elect the value of the spectroscopic mass of PG 2131+066 in- 
ferred by Miller Bertolami & Althaus (2006) for this comparison be- 
cause they use the same post-bom again PGl 159 evolutionary models 
we employ here in the determination of the asteroseismological mass, 
and because the spectroscopic masses quoted by Werner & Herwig 
(2006) are based on old post-AGB tracks. 



Quantity 



Spectroscopy Seismology This work 



Teff [kK] 
M, [Me] 

logg [cm/s^ 

log(L./Lo) 

log(R,/Re) 

Me„v [Mo] 

My [mag] 
Mbo\ [mag] 
Av [mag] 
d [pc] 
7T [mas] 



95 ±5 
0.58 + 0.1 
(0.55 + 0.1) 
7.5 + 0.5 



102.18+^;° 



0.608 ± 0.01 0.589+ 



1.41 + 0.5 
-1.73 

3.8 X 10-^^ 
7 69+"-^^ 

'•"^-0.18 



1.50 + 0.2 



7.63 _oi4 

1 c -7 +0.07 
^■-"-0.06 



+0.06 
-0.05 



-1.71 
0.016 



0.825^1^^ 
0.18 

830!^™ 

1 9+0.4 



Table 5. Same as Tableg] but for PG 1707+427. 



Quantity 



Spectroscopy Seismology This work 



[kK] 
M. [Mo] 

logg [cm/s-] 

log(L./Lo) 

log(R,/Re) 

Menv [Mq] 

Mv [mag] 
Mboi [mag] 
Av [mag] 
d[pc] 
n [mas] 



85 ± 4.5 
0.59 + 0.1 
(0.53 + 0.1) 
7.5 + 0.3 



0.57 ± 0.02 



1.36 



542+" '"'' 



7 CO +0.09 

1.40 + 0.04 
-1.68 + 0.04 
0.025 
7.25 ± 0.6 
1.25 ± 0.1 
0.12 

73Q+230 

1.4 + 0.4 



fit model has M, - 0.542 Mq, which is in agreement with the 
value derived from the average of the computed period spacings 
(M» ~ 0.566 Mq), but at odds (~ 9% lower) with that inferred 
from the asymptotic period spacing (M, = 0.597 Mq) (see Table 
[1). On the other hand, we note that M, for the best-fit model 
is in excellent agreement with the spectroscopic derivation of 
Miller Bertolami & Althaus (2006) (0.542 Mq versus 0.53 Mq), 
but is substantially lower than the spectroscopic value quoted by 
Werner & Herwig (2006) (0.59 Mq). Until now, the asteroseis- 
mological mass of PG 1707+427 has been more than 7% larger 
(AM» 0.04 Mq) than the spectroscopic mass if we adopt for the 
seismological mass the value found by Kawaler et al. (2004) and 
the derivation of Miller Bertolami & Althaus (2006) for the spec- 
troscopic mass. In light of our best-fit model, this discrepancy is 
strongly reduced to about 2% (AM, » 0.012 Mq). Finally, our 
best-fit model for PG 1707+427 is slightly more luminous than 
what is suggested by Kawaler et al. (2004). 
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4.3. The asteroseismological distance and parallax of 
PG 2131+066 and PG 1707+427 

As in our previous works — see, e.g., Corsico et al. (2007a) 
— we employ the luminosity of our best-fit models to infer the 
seismic distance to PG 2131+066 and PG 1707+427. Following 
Kawaler et al. (1995), we adopt BC = -6.0 + 0.5 for both stai's 
(Werner et al. 1991). We account for the interstellar absorption. 
Ay, using the interstellar extinction model of Chen et al. (1998). 
With all these ingredients the seismic distance, d, can be easily 
computed using the apparent magnitudes, which are my = 16.6 
for PG 2131+066 and = 16.7 for PG 1707+427 (Bond et 
al. 1984). We obtain a distance d ~ 830 pc and an interstellar 
extinction Ay ~ 0.18 for PG 2131+066 and d ~ 730 pc and 
Ay ~ 0.12 for PG 1707+427. 

Our estimation of the distance to PG 2131+066 is ~ 15% 
larger than that derived by CA06 (d ~ 116 pc). This is because 
our asteroseismological model is somewhat more luminous than 
thatofCA06 (log(L./Lo) = 1.57 versus log(L./LQ) = 1.37). 
On the other hand, our distance is 20 - 25% larger than that 
obtained by Reed et al. (2000) {d ~ 668 pc) on the basis of 
their own asteroseismological analysis. This difference can be 
understood on the basis that Reed et al. (2000) uses a lumi- 
nosity of log(L»/L0) ~ 1.4, somewhat lower than that of our 
best-fit model for PG 2131+066, of log(L,/LQ) ~ 1.6. On the 
other hand, our asteroseismological distance for PG 2131+066 
is about 1 .2 times longer than that quoted by Reed et al. (2000) 
of ~ 680 pc obtained on the basis of spectrum fitting, although 
both estimations are compatible at the Icr level. 

For PG 1707+427, our asteroseismological distance is in 
agreement with that quoted by CA06, of ~ 697 pc. Werner et al. 
(1991) obtain a distance to PG 1707+427 of- 1300 pc, substan- 
tially larger than our estimation, but still within the quoted error 
bars. The different value of Werner et al. (1991) is due to that 
they use a luminosity of log L^/Lq - 2.15, substantially higher 
than the luminosity of our best fit model (log L, /Lq - 1 .4). 

5. Comparison with the results of CA06 

Following the recommendations of an anonymous referee, we 
include in this section a detailed comparison between the 
PG1159 models and the asteroseismological results of the 
present paper and those of the previous study by CA06. 
These authors performed an asteroseismological analysis of four 
GW Vir stars (PG 0122+200, PG 1159-035, PG 2131+066, 
and PG 1707+427) on the basis of a set of twelve PG1159 
evolutionary sequences with different stellar masses (M* - 
0.53, 0.54, 0.55, ■ ■ ■ , 0.62, 0.63, 0.64M0) artificially derived from 
the full evolutionary sequence of O.5895M0 computed by 
Althaus et al. (2005). That sequence is one of the sequences 
we use in the present paper. Specifically, the sequences of 
CA06 were constructed using LPCODE by appropriately scal- 
ing the stellar mass of the O.589M0 sequence before the models 
reach the low-luminosity, high-gravity stage of the GW Vir do- 
main. Although this procedure leads to a series of unphysical 
stellar models for which the helium-burning luminosity is not 
consistent with the thermo-mechanical structure, the transitory 
stage vanishes shortly before the star reaches the evolutionary 
"knee" in the HR diagram (see Fig. 2 of CA06). As a conse- 
quence, those PG1159 models were not suitable for the high- 
luminosity, low-gravity regime corresponding, for instance, to 
RX J21 17.1+3412, NGC 1501, K 1-16, HE 1429-1209, etc. 

Because the sequences of CA06 with different stellar masses 
were created starting from a single sequence with M, - 
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Fig. 4. Comparison between some PG1159 evolutionary tracks 
employed in the present work (thick lines) and those of CA06 
(thin lines). The values are in solar masses. Note that the track 
corresponding to the sequence of « O.59M0 is the same for the 
two sets of computations. 



O.589M0, the central abundances of C and O and the spatial 
extension of the C-O core are not completely consistent with 
the value of the stellar mass, except in the case of the models 
of the O.589M0 sequence itself. For instance, a model of CA06 
with M, - O.64M0 has a C-O core that is somewhat smaller 
and the central abundance of O is substantially higher than what 
would be expected if the complete evolution of the progeni- 
tor star were performed, as it is the case for the models em- 
ployed in the present paper. Fig. 2 of Miller Bertolami & Althaus 
(2006) clearly illustrate that, when the complete evolution of the 
PG1159 progenitor stars is taken into account, different stellar 
masses are associated with different central abundances of C and 
O, and different sizes of the C-O core. Specifically, the more 
massive the models, the lower (higher) the central abundance of 
O (C) and the larger the C-O core. In summary, for a given value 
of M», and /or stages after the evolutionary knee, the only struc- 
tural/physical difference between the PGl 159 models employed 
in the present work and those of CA06 is related to the size of 
the C-O core and the central abundances of O and C. 

In Fig. |4] we compare some evolutionary tracks (M, = 
0.530, 0.542, 0.589, 0.609Mo) of the PGl 159 sequences em- 
ployed in the present work with the corresponding evolutionary 
tracks of CA06 (M» = 0.53, 0.54, 0.59, 0.6IM0). A careful in- 
spection of this figure reveals that both sets of tracks generally 
differ, but when the models reach the beginning of their white 
dwarf stage (log Teff < 5), they turn be in very close agreement. 
This agreement is reached earlier in the case of sequences with 
stellar masses close to M» = O.589M0, the value of the sequence 
from which the remainder sequences of CA06 were generated. 
Interestingly enough, the regime in which the tracks of CA06 are 
in agreement with those employed in the present paper embraces 
the location of PG 2131+066 and PG 1707+427, which have 
log Teft - 4.97 and log Tetf - 4.93, respectively. Because of this, 
it is expected that the global pulsation properties (i.e., asymp- 
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Table 6. Comparison between the stellar mass values (in M©) 
obtained in CA06 and in the present work. 



r\j 21 jl+UOO 




A TT 


Period fit 


CA06 


0.615 


0.575 


0.60 


This work 


0.627 


0.578 


0.589 


AM./M. 


2 % 


0.3 % 


2 % 


PG 1707+427 


Anj 


att" 


Period fit 


CA06 


0.595 


0.565 


0.55 


This work 


0.597 


0.566 


0.542 


AM./M. 


0.3 % 


0.1 % 


1.5 % 



totic period spacing, average of the computed period spacing) of 
both sets of models in that regime should nearly agree, and con- 
sequently the asteroseismological inferences on PG 2131+066 
and PG 1707+427 based on these two different sets of models 
should not be substantially distinct. 

In Table |6] we present a comparison between the stellar 
masses inferred in CA06 (their Table 1) and in the present study. 
Note that, not surprisingly, the stellar masses derived in CA06 
are in excellent agreement with the values obtained in the present 
work, irrespective of the particular method employed, being the 
differences in all of the cases below 2%. This is a clear indica- 
tion that the asteroseismological results of this paper and that of 
CA06 for PG 2131+066 and PG 1707+427 are not seriously af- 
fected by the differences between both sets of models. This adds 
credibility and robustness to the asteroseismological results of 
the present study. This conclusion should change for the case 
of PG1159 pulsators with stellar masses too departed from the 
value ~ O.589M0 and/or located at earlier evolutionary stages. 

6. Summary and conclusions 

In this paper we presented an asteroseismological study 
of the high-gravity, low-luminosity pulsating PG1159 stars 
PG 2131+066 and PG 1707+427 and of the high-luminosity 
PNNV [WCE] star NGC 1501. This is the fourth article of a 
short series of studies aimed at exploring the internal structure 
and evolutionary status of pulsating PG1159 stars which have 
been intensively observed through multi-site campaigns. Our 
analysis is based on the full PG1159 evolutionary models of 
Althaus et al. (2005), Miller Bertolami & Althaus (2006) and 
Corsico et al. (2006). These models represent a solid basis to 
analyze the evolutionary and pulsational properties of pre-white 
dwarf stars like PG 2131+066, PG 1707+427, and NGC 1501. 

We first used the observed period-spacing data to obtain esti- 
mations of the stellar mass of PG 2131+066, PG 1707+427, and 
NGC 1501. The results are summarized in Table [T] where we 
provide a summary of results from the present and the previous 
works by us, and also from other pulsation and spectroscopic 
studies. We obtained three mass values for each star: the first 
one by comparing the observed period spacing with the asymp- 
totic period spacing of our models (an inexpensive and widely 
used approach that does not involve pulsational calculations); 
the second one by comparing the observed period spacings with 
the average of the computed period spacing (an approach that re- 
quires of detailed period computations); and the third one on the 
basis of the approximate formula of Kawaler & Bradley (1994), 
which is based on the behavior of the asymptotic period spac- 
ing of a large grid of quasi-evolutionary PG1159 models. The 
first and the third approaches are almost equivalent, and lead to 



similar, somewhat overestimated values of the stellar mass for 
PG 2131+066 and PG 1707+427. The second approach, clearly 
more realistic, conducts to smaller values of M» (and closer to 
the spectroscopic inferences) in the case of PG 2131+066 and 
PG 1707+427, and virtually the same M. value than that ob- 
tained from the first approach in the case of NGC 1501. 

In the second part of our work, we sought for the mod- 
els that best reproduce the individual observed periods of each 
star. The period fits were made on a grid of PG1159 models 
with a quite fine resolution in effective temperature (AT^ff ~ 
10 - 30 K) although admittedly coarse in stellar mass (AM» ~ 
0.01 - 0.08 Mq). We found asteroseismological models only for 
PG 2131+066 and PG 1707+427. For NGC 1501 we were un- 
able to find a clear and unambiguous seismological solution 
due to the existence of numerous and equivalent minima char- 
acterizing the quality function employed in the period-fit pro- 
cedure. The pulsational properties of the "best-fit" models for 
PG 2131+066 and PG 1707+427 are summarized in Tables |2] 
and [3] respectively. In particular, we predict the values of the 
rates of period change to be positive and in the range (2 - 5) x 
10"'^ s/s. Unfortunately, we have been unable to check the real- 
ity of this prediction because the lack of any measurement of n 
for PG 2131+066 and PG 1707+427 for the moment. The struc- 
tural characteristics of these best-fit models are shown in Tables 
|4]and|5] In particular, the seismological masses are closer to the 
spectroscopic ones in light of our best-fit models, as we found in 
our previous works. In these tables we also show the seismologi- 
cal distances and parallaxes of PG 2131+066 and PG 1707+427. 
We found a reasonable agreement between our results and those 
of Kawaler et al. (1995), Reed et al. (2000), Kawaler et al. 
(2004), and CA06. We stress that almost all differences between 
our results (Sections 14.21 and l43b and those of earlier works are 
within the quoted errors. 

In summary, in this work we have been able to estimate the 
stellar mass of PG 2131+066, PG 1707+427, and NGC 1501 
on the basis of the period-spacing information alone. We have 
also been successful in finding asteroseismological models for 
PG 2131+066 and PG 1707+427 from period-to-period com- 
parisons. In particular, the T^ff and \ogg of the best-fit mod- 
els are in very good agreement with the spectroscopic mea- 
surements. Unfortunately, we fail to found an asteroseismolog- 
ical model for NGC 1501. In principle, this shortcoming of our 
study could be indicating some inadequacy inherent to the stel- 
lar modeling. Another possible alternative could be the fact that 
the period spectrum of NGC 1501 includes periodicities asso- 
ciated with g-modes with radial orders very spaced from each 
other, in such a way that our procedure of period-fit is inef- 
ficient to isolate a clear and unambiguous asteroseismological 
solution. On the other hand, it would be kept in mind that while 
both PG 2131+066 and PG 1707+427 are classified as PG1159 
stars, NGC 1501 is a [WCE] star Although both classes are sus- 
pected to form an evolutionary sequence, this possibility is still 
under debate (Crowther 2008; Todt et al. 2008) and it could not 
be the case. Therefore, the failure of our models to fit the pe- 
riod spectrum of NGC 1501 might be indicative that this star 
has a very different evolutionary history than PG 2131+066 and 
PG 1707+427. 

We have also included a comparison between the models and 
results of the present work and those of the study by CA06. The 
models employed in the present work are the result of the com- 
plete evolution of the progenitor stars, and as a result, they are 
characterized by central chemical abundances and a size of the 
C-O core which are consistent with the value of the stellar mass. 
This is not the case for the models employed by CA06, which 
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were artificially derived from the full evolutionary sequence of 
0.589Mo computed by Althaus et al. (2005). In spite of these dif- 
ferences, our asteroseismological results are in excellent agree- 
ment with those of CA06. This adds credibility and robustness 
to the results of the present study. 

As the main conclusions of the present work, we can men- 
tion: 

- The full evolutionary models of PG1159 stars employed 
in the present work lead to asteroseismological results for 
PG 2131+066 and PG 1707+427 that do not diff^er substan- 
tially from those predicted by CA06 on the basis of evolu- 
tionary sequences generated artificially. We note, however, 
that this agreement between both sets of computations is 
valid for PG1159 stars located at the low-luminosity, high- 
gravity regime after the stars have passed the evolutionary 
knee in the HR diagram. It should be kept in mind that this 
conclusion should change for the case of stars located at ear- 
lier stages of evolution. 

- At present, the PG1159 evolutionary models used in this 
work — and in the previous studies of this series — remain 
the only suitable for asteroseismological inferences on stars 
that are located at the high-luminosity, low-gravity regime 
before the evolutionary knee, such as RX J21 17.1+3412 and 
NGC 1501. 

- The detailed fitting of the individual periods (Sect. gives 
somewhat different masses than analysis based on asymp- 
totic period spacing (methods 1 and 3 of Sect. [3]), but in very 
good agreement with the values of M, derived from the av- 
erage of the computed period spacings (method 2 of Sect. 

0. Thus, method 2 is a very appropriate way to estimating 
stellar masses, and detailed period fits do not significantly 
improve the mass determinations. We note, however, that 
the period-fit approach yields an asteroseismological model 
from which one can infer, in addition to M», the luminosity, 
radius, gravity, and distance of the target star In addition, the 
period-fit approach does not require — in principle — exter- 
nal constraints such as the spectroscopic values of T^s and g, 

1. e., the method works "by letting the pulsation modes speak 
for themselves" (see Metcalfe 2005 for an interesting discus- 
sion about this). 

- The nonadiabatic stability analysis does not at the moment 
predict instability for all of the fitted modes. This means 
that, in the frame of the linear nonadiabatic pulsation the- 
ory, some pulsation modes detected in PG 2131+066 and 
PG 1707+427 should not be excited. It is not clear at this 
stage the origin of this discrepancy. Maybe it could be at- 
tributed to the extreme sensitivity of the stability analysis of 
PGl 159 stars to the exact amounts of the main atmospheric 
constituents (see Quirion et al. 2004 for details). 

- The main conclusion of this series of papers — 
the present work and Corsico et al. (2007a, 2007b, 
2008) — is that for most well-observed pulsating PGl 159 
stars (RXJ21 17.1+3412, PG 0122+200, PG 1159-035, 
PG 2131+066, and PG 1707+427) it is possible to found a 
stellar model (the asteroseismological model) with M» and 
Teff near the spectroscopic measurements to a high internal 
accuracy. The next step is of course an assessment of the 
question if the asteroseismological models can provide more 
accurate masses for these objects. The scatter in the masses 
derived from the different asteroseismological methods (see 
TablelTJ suggests that it may not be the case. In fact, when all 
asteroseismological methods are considered, the uncertainty 
in the determination of the mass amounts to ~ O.O5M0, com- 



parable to the spectroscopic one (~ 0.05 - O.IM©; Werner 
et al. 2008). However, it is worth noting that, when results 
based on asymptotic period spacing (an approach that is 
not correct for the high-gravity regime of PGl 159 stars; see 
Althaus et al. 2008a) are not taken into account, the scatter- 
ing in the derived masses is of only ~ O.O2M0. 

We close the paper by noting that the PGl 159 evolutionary 
models employed in our series of asteroseismological studies 
are characterized by thick He-rich outer envelopes, as they are 
predicted by the standard theory for the formation of PGl 159 
stars. However, Althaus et al. (2008b) have recently demon- 
strated that the assumption of thinner He-rich envelopes solves 
the longstanding discrepancy between the measured rates of pe- 
riod change in the prototypical star PG 1159-035 and the pre- 
dictions of theoretical models. In view of this important re- 
sult, we are planning future asteroseismological studies for all 
the pulsating PGl 159 stars analyzed in our series of articles, 
but with non-canonical PGl 159 models characterized by thin- 
ner He-rich envelopes. 
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